Phonon renormalization from local and transitive electron-lattice couplings in 

strongly correlated systems 



On 
O 
O 
(N 
•*-> 
O 

O 

(N 



E. von Oclscn, 1 A. Di Ciolo, 2 J. Lorenzana, 2 ' 3 G. Seibold, 1 and M. Grilli 2 

'Institut fur Physik, BTU Cottbus, PBox 101344, 03013 Cottbus, Germany. 
2 SMC-INFM-CNR, Dipartimento di Fisica, Universita di Roma "La Sapienza", P. Aldo Mow 2, 00185 Roma, Italy. 

3 ISC-CNR, Via dei Taunni 19, 00185 Roma 
(Dated: October 30, 2009) 

Within the time-dependent Gutzwiller approximation (TDGA) applied to Holstein- and SSH- 
Hubbard models we study the influence of electron correlations on the phonon self-energy. For 
the local Holstein coupling we find that the phonon frequency renormalization gets weakened upon 
increasing the onsite interaction U for all momenta. In contrast, correlations can enhance the phonon 
frequency shift for small wave-vectors in the SSH-Hubbard model. Moreover the TDGA applied to 
the latter model provides a mechanism which leads to phonon frequency corrections at intermediate 
momenta due to the coupling with double occupancy fluctuations. Both models display a shift of the 
nesting-induced to a q — instability when the onsite interaction becomes sufficiently strong and 
thus establishing phase separation as a generic phenomenon of strongly correlated electron-phonon 
coupled systems. 

PACS numbers: 71.10.Fd,71.38.-k, 74.72.-h 
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I. INTRODUCTION 

Transition metal compounds are usually character- 
ized by strong electron-electron and electron-phonon 
interactions (for an overview cf. Ref. [l]). The interplay 
of these interactions can give rise to a large variety of 
interesting electronic properties which also reflect in 
the energy and momentum structure of the phonons. 
In this regard a well known phenomenon is the oc- 
curence of bond-stretching phonon anomalies which 
are observed in high-T c cuprates La^^Sr^CuO^^ 
and YBa 2 2Cu 3 6+x ^ HgBa 2 Cu0 4+ , 5 r r - 

Ndi. 86 Ceo.i4Cu0 4+ 5, 8 Bi 2 Sri. 6 La .4Cu 2 O 6 +,?2 but 
also in Pb- and K-doped BaBi0 3f ^ Sr 2 Ru0 4 ^- 
Lai.69Sr .3iNiG\fi2 and Lai-^Sr^MnOaii 3 - 

The interpretation of such experiments requires an 
understanding of the renormalization of phonons in a 
strongly correlated electron system which is the purpose 
of the present paper. Such an analyis is usually based on 
Hubbard-type models where the two most popular vari- 
ants to implement an electron-lattice coupling arc cither 
to consider the dependency of onsite energies or of the 
hopping integrals as a function of some atomic coordi- 
nates. In the first case, when one restricts to the in- 
teraction between electron density and coordinate of the 
same lattice site, the coupling is usually termed a Hol- 
stein or molecular crystal model^ In the second case, 
when the hopping between nearest-neighbor sites is ex- 
panded in terms of the positions of these sites the re- 
sulting electron-lattice coupling is often named after Su, 
Schrieffcr and Heeger (SSH) which have used this type of 
interaction for the analysis of solitons in polyacetylene:- 5 

Since lattice vibrations are partially screened by the 
electrons the latter can have a profound influence on 
the effective dispersion of the phonons. One example 
is the Kohn anomaly^ caused by the abrupt change of 
electronic screening near wave- vectors q which are twice 



the Fermi momentum fc^ (nesting). Moreover, electronic 
correlations alter the electron dynamics of the system 
and thus also affect the phonon dispersion. In case of a 
Holstein coupling, and due to the fact that correlations 
reduce the charge correlations at kp in favor of an en- 
hancement of the q = response, it has been show n 17 i 18 
that the formation of a charge density wave (CDW) is 
suppressed in favor of a phase separation instability. For 
the phonons this has the immediate consequence that the 
softening is shifted from the Kohn anomaly wave-vector 
q = 2kp to q = 0. 

Naturally the interplay between lattice and elec- 
tronic degrees of freedom has already been investi- 
gated by means of several techniques. Some of the ap- 
proaches, like quantum Monte CarlO ) 19 ' 20 i 21 ' 22 i 23 ' 24 i 25 ex- 
act diagonalization ) 26 i 27 i 28 i 29 and Dynamical Mean Field 
Theory (DMFT)2iE^ili2i^ are intrinsically non- 
perturbative but are numerically challenging and suffer 
some limitations (small lattice sizes, no momentum res- 
olution, etc.). On the other hand (semi) analytical ap- 
proaches like variational ones, slave boson and large- 
N expansionil^ 3 ^^^^ 3 -^^ deal with infinite 
systems, but within approximate treatments. It should 
be noted that most of the work has been done on the 
Hubbard-Holstein model. The papers which explicitely 
deal with an Hubbard-SSH mode&^L^dL^. mainly 
focus on the interplay between correlations and dimcr- 
ization or superconductivity, respectively. 

In this paper we want to study the correlation effects 
on the phonon excitations for both Holstein and SSH 
models on the same footing. To this aim we need a 
method which is not numerically very demanding, but 
still provides a quantitatively acceptable treatment of 
the strongly correlated regime. In this regard we find 
the Gutzwiller approximation (GA) supplemented with 
RPA-type fluctuations, the so-called time-dependent GA 
(TDGA), a good compromise. This technique corre- 
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sponds to Vollhardt's Fermi liquid approach^ 7 - with the 
fluctuations extended to finite frequencies and momenta. 
Electron-phonon interactions will be treated in the Born- 
Oppenhcimer approximation. 

It is worth mentioning that the TDGA approach has 
been tested in various situations and found to be accurate 
compared with exact diagonalizationi 48 i 49 ' 50 i 51 Compu- 
tations for realistic models have provided a description 
of different physical quantities in agreement with exper- 
imentW>£* 

In the discussion of results we will restrict to one- 
dimensional systems. This clearly simplifies the compu- 
tations and the presentation of the results. It has the 
drawback that strictly one dimensional systems are those 
for which our Fermi liquid like approach is expected to 
be less suited. Thus our results should not be taken too 
literally in this case. On the other hand the qualitative 
behavior we find is rather independent of dimension. For 
example for the Holstein model we find here the same 
qualitative behavior for the charge response as we have 
found before in larger dimensions. 63 Despite the inade- 
quacy of a Fermi liquid treatment for strictly one dimen- 
sional systems several real systems are only quasi one 
dimensional and in those cases our computations can be 
applicable. We also found that in certain filling ranges 
the results are surprisingly accurate. 

The scheme of our paper is as follows. In Sec. II we 
define the model and show how the Hubbard-Holstein 
and -SSH hamiltonians are represented within the GA. 
After introducing the TDGA in Sec. IIC the phonon self- 
energy for both Holstein- and SSH-couplings is derived in 
Sec. IID,E. In Sec. IIIA we first investigate the quality of 
our technique by comparing with exact and Monte-Carlo 
results in one dimension and furtheron present the results 
for the phonon self-energies of the Holstein- and SSH- 
Hubbard model in Sees. IIIB,C. Our conclusions can be 
found in Sect. VII and the details of our calculations are 
given in the Appendix. 



II. FORMALISM 
A. Model 

Our investigations are based on the following hamilto- 
nian 

H = H e + H e - P h + H p h (1) 

where H e denotes the Hubbard model, H e _ p h the cou- 
pling between electrons and phonons and H p h the bare 
phonon part. Here we restrict to one-dimensional sys- 
tems and consider for the electronic part hopping be- 
tween nearest neighbors 

i,cr i 

(2) 



where c\ l destroys (creates) an electron on lattice site 

Ri and n iy(7 = cj ;0 .c ii0 -. 

We consider two types of electron-phonon coupling. 
The first is a local Holstein interaction initially motivated 
from a molecular crystal model 

n™ h = -aJ2n(n i ,„-(n i , ( ,)) (3) 

i,a 

where Ui is the coordinate of an internal mode of the 
molecules affecting the site energy at site i. Its dynamics 
is described by 

i i 

and corresponds to dispersionless lattice modes. Here K 7 
M denote elastic constant and reduced mass and pi are 
the conjugate momenta at site i. 

The second type of coupling originates from the de- 
pendence of the electronic hopping on the atomic coor- 
dinates. For the one-dimensional model under consid- 
eration this so-called SSH or Peierls interaction is given 
by 

i,(T 

(5) 

and in this case the lattice dynamics is determined from 

H Ph = l K J2( u ^ - + ^T,pi ( 6 ) 



B. Gutzwiller approximation 

We treat the model Eq. {T]) within the Gutzwiller 
approximation (GA) supplemented by Gaussian fluctua- 
tions the evaluation of which are outlined in the next sec- 
tion. The GA can either be motivated from a slave-boson 
approach^ or a variational Ansatz formally evaluated in 
infinite dimensions?^ The variational wave function is 
given by \^f) = P\<j>), where the Gutzwiller projector P 
acts on the Slater determinant \<fi). This approach in- 
corporates the correlation induced renormalization of the 
kinetic energy and treats the Hubbard on-site interaction 
via the variational double occupancy parameters Dj. For 
each term in the Hamiltonian we derive an energy func- 
tional E GA [p,D] = (^\H\^f) where we have introduced 
the one body density matrix associated with the Slater 
determinant p% t j,cr = {4>\ c \ a c j,a\4>)- Specifically the en- 
ergy functional for the electronic part Eq. @ reads as 

E° A [p, D] = -t Z,>Zi+l, CT (Pi,i+l,a- + Pi+l,i,a) 

+ UY J D 1 (7) 
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and the z-factors are given by 



z _ \J (Pit.a - A)(l - Pit + A) + \/{pa-<j ~ Di)D i 

\J Pii,a{X — Pii. a) 

(8) 

Further on p^ = J2 a Py>- 

Within the Born-Oppenheimcr approximation the 
electronic expectation value of H e ^ p h determines the lat- 
tice potential. In contrast to the local Holstein coupling, 
for electronic degrees the transitive SSH electron-phonon 
interaction Eq. is also renormalized by the z-factors 
and the corresponding energy functional reads as 



E 



ssh,GA 
e—ph 



14+1 - Ui) z it<T z i+ i t , 



x (Pi,i+1,<7 + Pi+l,i,a) ■ 



(9) 



The GA variational ground state is then obtained upon 
minimizing E GA with respect to {D}, {u}, and p under 
the constraint that the latter derives from a Slater deter- 
minant, i.e. p 2 = p. In the following our starting point 
will be an homogeneous state. The formation of possible 
charge-density wave (in case of the Holstein model) and 
dimerized (in case of the SSH model) states will appear 
as instabilities of the homogeneous state. Thus our ini- 
tial ground state is characterized by 14 = 0, pu. a = po/2, 
Di = Dq and Zj = zq. The ground state energy per site 
is therefore simply determined from Eq. ([7]) and reads 



E GA [p,D]/N = z 2 e Q + UD 



(10) 



and eo denotes the energy per site of a non-interacting 
system with charge density pg. For later use we also 
denote the critical value of the onsite repulsion U = U c — 
32t/7r for a half-filled one-dimensional chain at which the 
Brinkman-Rice transition (i.e. complete localization of 
the charge carriers with Dq — 0) takes place. 



C. Time-dependent Gutzwiller approximation 

In order to study fluctuations beyond the GA saddle- 
point solution, necessary for the evaluation of the phonon 
self-energies, we use the time-dependent GA (TDGA) 
which has been developed in Refs. |48||49| . 

We briefly illustrate the formalism for the electronic 
part H e and furtheron show how lattice fluctuations can 
be implemented into the theory. Further details can be 
found in Refs. [HlHil 

We study the response of the system to a small time 
dependent external field which produces time dependent 
fluctuations in the density matrix Sp and the double oc- 
cupancy SD. This can then be obtained by expanding 
E GA up to quadratic order in the density and double 
occupancy fluctuations Sp and SD. 

For a translationally invariant ground state it is con- 
venient to perform the expansion in momentum space. 



Besides the local fluctuation Spq, we introduce the bond 
charge fluctuation^, 

STi =~t )^ (Spi+,,,i,cr + Spi,i+ V ,<j)- 

It is convenient to introduce the hopping factor in the 
definition of STi so that it can also be interpreted as a 
local kinetic energy. Notice however that the z factors 
are omitted. The Fourier transform is given by 

ST q = -2i^[cos(fc + c/) +co8(k)}Sp k+q>kt<7 . (11) 

k,a 

For latter use it is also convenient to introduce the anti- 
symmetric combination of bond charge fluctuations: 

6Tf = -t 1](5p i+ , h i,cr + Sp hi+V ^). 

ar\— ±1 

with Fourier transform: 



5T~ = -2it^T [siu(fc + q)- sin(ft)] Sp 



k-\-q.k,a • 



(12) 



k.a 



It is easy to check that the two fluctuations are related 
by a function of q: 



st: 



i tan 



(1) 



5T a 



(13) 



The second order-energy expansion in the charge channel 
of Eq. follows as: 



e: 



CM, (2) _ 



N 



\ V q S Pq$P-q + Z ^ Z D SD q ST^ 



+ 2 Z °( Z> + Z +-)^2 S Pq ST -q 

q 

+ ]T L q Sp q 5D_ q + 1^2 U q SD q 5D^ 



(14) 



with the following definitions: 



V q 
L„ 



e z 



{z'l + + 2z'l_ + z"__) + -{z' + z' + _) 2 e cos(g) 



eoZo{z+ D + z"_ D ) + z' D (z' + z' + _)e cos(q) 



U q = 2e z z'o + 2 ( z 'd) e ocos(<?) 

where z' and z" denote derivatives of the hopping factors 
which are given in the Appendix. 

The double occupancy fluctuations can be expressed 
in terms of the density fluctuations by use of the an- 
tiadiabaticity condition which assumes that the double 
occupancy adjust to the instantaneous configuration of 
the charge^ 



dE, 



CM, (2) 



dSD„ 







(15) 
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T 



I r 

- n=l 

— n = 0.98 

- ■ n = 0.9 
■ - n = 0.8 



where 



U[U] 



FIG. 1: (Color online) The interaction coefficient A q= o for 
various values of particle density n in units of the Brinkman- 
Rice critical onsite repulsion(7 c for a one-dimensional chain. 



A q B q 

B* C q 



is the interaction kernel. The elements of W q e e 



(17) 



are given 



by 



T 2 

A — V — _2 

Ug 



(18) 



B„ 



C q = 



1 



z ( z' + z' + _ - z' L 



2L, 
U„ 



4i^ C os 2 («/2) 



and the long- wave limit of the (11) element A q=0 , which 
dominates the interaction kernel close to half-filling is 
shown in Fig. Q] For small U one finds A q= o w £7/2 
whereas A q= o is enhanced close to the Brinkmann-Rice 
transition U = U c . At exactly half- filling one finds 



U(U C + U){U -2U C ) 
4U C (U - U c ) 



(19) 



and one obtains the following functional which only de- 
pends on the local and intersite charge deviations: 



E GA,(2) = J_ V 



6T 



A q B q 
B*q C q 



ST- 



(16) 



for U < U c . 

Since the energy expansion in Eq. (116j) is a quadratic 
form in 5p q and 5T q (see also Eq. (ITTj) ). it is useful to 
introduce the following susceptibility matrix for the non- 
interacting system \q : 



ka 



-2t [cos(fe) + cos(fc + q)] 



-It [cos(fc) + cos(fc + q)] At [cos(A;) + cos(fc + q)] J ui + eu+ q - ek — iO 



nk+q,a — Tike 



(20) 



r 



The susceptibility for the interacting system is then ob- 
tained from the following RPA series 



Xq = Xq + 



Xq°Wq £ 



'Xq 



(21) 



where the element (x g )n corresponds to the correlation 
function for the local charge response. 



D. Phonon self-energy for the Holstein coupling 

As mentioned above, due the local nature of the Hol- 
stein coupling Eq. ([3]) it is not renormalizcd by the z- 
factors and thus its quadratic contribution to the energy 
expansion is given by 



E 



hoi, (2) 
e—ph 



1 



\ " 



&Pq 



(22) 



where Q q = J^i exp(—iqRi)iii denotes the Fourier trans- 
formed (normal) coordinate fluctuation (remember that 
the saddle-point solution has Uj = uq = so that we can 
skip the 5 symbol). 

Similar we write the lattice energy Eq. (H|) as 



E 



hol,(2) 
ph 



-E 

2N ^ 



P q P-q 

M 



MV?Q„Q_ 



(23) 



with n 2 = K/M. 

The small time dependent deviation from the elec- 
tronic ground state 8p q will act as a force on the lattice 

coordinates via -E^-ph an< ^ the corresponding equation 
of motion reads 



MQ q + MtfQq 



-N- 



dE. 



hoi, (2) 
e—ph 



dQ- 



aSpq 



(24) 
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As a consequence the lattice vibrations are shifted to 
new frequencies oj q which depend on the electronic charge 
fluctuation 5p q 



n 2 ) Q q 



a X 
M 



(25) 



On the other hand 5p q can be determined from linear 
response theory when we view Eq. (|22|) as a small per- 
turbation on the electronic system. With the charge sus- 
ceptibility derived in the previous section one has 



s Pq = (Xq)u(-®Qq) 
which upon inserting in Eq. (|25p yields 



lo 2 = n 2 



M 



(x. 



q)ll 



n 2 



2fffi„ 



(26) 



(27) 



where T, q denotes the phonon self-energy. In these equa- 
tions Xg( w ) an d ^g(w) should be evaluated at lu = uj q . 
Since the phonon dynamics is assumed to be much slower 
than the electron dynamics it is a good approximation to 
evaluate the susceptibilities in the static limit. 

We see that in the case of the Holstein coupling the 
phonon frequency shift is solely determined by the local 
charge susceptibility renormalized by electronic correla- 
tions within the GA. For later use we define the ratio 
between self-energies in the correlated and uncorrelated 
case 



r„ 



£ 9 (*7 = 0) 



which for the Holstein coupling becomes 



r 9 = 



(x q ) 



q)ll 



(x2) 



ii 



(28) 



(29) 



In addition it is convenient to introduce the following 
coupling constant which has energy units 



9hol = a 



1 



2MQ 



(30) 



so that the phonon self-energy for the Holstein model is 
given by 



S 9 — 9hol(Xq)ll- 



(31) 



E. Phonon self-energy for the SSH coupling 

The case of the SSH electron-lattice coupling is more 
subtle since the interaction energy Eq. ([9]) depends on 
the z-factors the fluctuations of which we have to consider 
in the evaluation of T, q . 

Expansion of Eq. §§§ up to second order in the fluc- 
tuating fields yields for the Fourier transformed effective 



electron-lattice interaction 



E 



ssh,(2) 
e—ph 



aZ oiy"j T q Q-q 



q 

+ iae Q Zo(z' + z +-)-^ sm(q)6p q Q- q 
q 

+ 2iae z z' D — ^ sm(q)SD q Q- q 



(32) 



with the same definitions already introduced in Sec. Ill CI 
Besides the coupling to the transitive fluctuations 5T q 
the correlations induce a coupling of the lattice to local 
density (5p q ) and double occupancy (SD q ) fluctuations. 
Unlike the Holstein case the antiadiabaticity condition 
now includes the electron-lattice coupling Eq. (|32|) in 
addition to the bare electronic part Eq. JI| 



d 



E 



GA,(2) 



E 



sh,(2) 



ph 



d5D a 



= 



(33) 



so that the double occupancy fluctuations can be ex- 
pressed via the density and lattice fluctuations: 



dD q = 2iaz z D — F} — Q q - j^-dp q tz^oT,, 



U a 



U 



U a 



(34) 



Inserting Eq. (f3"4")l into Eqs. (|14l32p and including also 
the (fourier transformed) lattice part Eq. ^ yields 



E 



ssh,(2) _ E GA,{2) 



E 



sah,{2) 



E ssh,(2) (35) 



tot -^e ~ '-'e—ph ~ ph 

where E^ A ' {2) was derived in Sec[HCl Eq. (JTSJ). The 
effective coupling of the lattice to the electronic density 
fluctuations is given by 

E t ' (2) = a \ E K*P*Q-* + <4 E W 2 5T q Q_ q 



N 



where we used Eq. (I13[) to eliminate the antisymmetric 
fluctuations and introduced the elements of the vector 

■yy^el-ph 



W 1 



ie z sin(g) 



z' + z' 



2z'J^ 

D u 



(37) 



= iz£ttm{q/2)-2ie Q (z z' D ) 



, N 2 sin (<?) 



U„ 



The lattice part becomes 



E 



h,(2) _ J_ [ P q P -q 

ph 2N ^ \ 2M 

q v 



Mn 2 q Q q Q- q } . (38) 



Interestingly the elimination of the double occupancy in- 
troduces a novel renormalization of the phonon disper- 
sion 



K 



1 



, ,2 sin 2 (g) 



n 2 q = 2^ (1 - cos(g)) - — (2ae z Q z' D y (39) 
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The (squared) dispersion is composed of the acoustic 
branch ~ 1 — cos(q) of the uncorrelated atomic chain 
and a contribution which arises from the double occu- 
pancy fluctuations. To the best of our knowledge this 
is the first time such correlated renormalization of the 
phonons is reported. We will show below that it induces 
a phonon softening in addition to the contribution which 
comes from electronic screening. 

The phonon self-energy in the present case can be de- 
rived via a similar procedure as before, taking into ac- 
count the vectorial character of the susceptibility and 
reads 



S 9 = 9a 



w q e '~ p/l 



Xq 



w, 



el—ph 



(40) 



where the electronic susceptibility matrix \ q is obtained 

from Eq. (f2"Tj) . We also defined the dimensionless cou- 
pling constant 



9q 




2M 



(41) 



(42) 



In the uncorrelated limit U = the self-energy reduces 



to 



S 9 (C/ = 0) = 4 5? 2 ([/ = 0)sin 2 (g/2)( X «) 



22 



(43) 



where the coupling constant Eq. (|41"j) has to be eval- 
uated with the bare acoustic dispersion fl q (U = 0) = 
2* (l-cos(g)). 

Finally it should be noted that the parameters of the 
problem are given by hojQ = yj K/M , g and U measured 
in units of the hopping t = 1 and we set K = 1 in the 
following. 



where vp = v F z 2 is the quasiparticle Fermi velocity 
and v p = Vpy/l + 4:A /(ttvf) is the velocity of the 
(quasi)particle-hole excitations with A = A q=Q defined 
in Eq. (fTS]) . The compressibility k = — (x g ^o)n( w = 0) 
follows as 



2v F 

TTV 2 , 



1 



IT VF 



■ 4^0/71- 



(45) 



In the weak coupling limit these expression coincide with 
the perturbative expressions for the Tomonoga-Luttingcr 
liquid. 58 As in the exact case the effective interaction and 
the Fermi velocity gets renormalized upon increasing U. 




2 


\ 


b 
\ 


■ — u/u c =o 


1.5 




U/U c =0.2 




£ 1 




\ 


U/U c = l 


0.5 






U/U c =oo 











0.2 



0.4 0.6 
n 



0.8 



III. RESULTS 

A. Properties of the TDGA and comparison with 
exact results 

Since in this paper we restrict ourselves to one- 
dimensional systems it is necessary to analyze the quality 
of our approximate TDGA scheme in this special limit. 
With regard to the following investigation of phonon self- 
energies it is important to figure out which features can 
be considered as generic and also show up in higher di- 
mensions where our mean-field type approach is expected 
to perform better. 

Since the phonon self-energies are completetly deter- 
mined by the charge susceptibility (at least in the Hol- 
stein case) (xq)u we first analyze the corresponding 
TDGA result which for small wave- vectors q and close 
to half-filling can be expanded as 

(X 9 )n « 2 - 2 VF f ]2 (44) 
7T lu 2 - (v p q) 2 



FIG. 2: (Color online) Id-Charge compressibility as a func- 
tion of n and for different values of U/U c calculated with the 
Bethe ansatz - exact Id solution (a) and with the TDGA 
(b). Here U c ~ 32t/ir is the Coulomb repulsion at which the 
Brinkman-Rice transition takes place for n = 1 in the GA. 
In the exact solution the metal insulator transition occurs at 
U — therefore U c = 32t/ir is used only as an energy unit. 

In Fig. [2] we compare the TDGA charge compressibility 
with the exact Id solution of the Hubbard mode h 59 ' 60 i 61 
The renormalization of vf and the effective interaction 
pushes the qualitative agreement with exact results to 
larger values of U than the traditional HF+RPA ap- 
proach. Strong differences arise close to n = 1. As soon 
as the interaction is switched on the exact compressibility 
diverges. This can be understood in the strong coupling 
limit where the charge degrees of freedom can be mapped 
to a spinless Fermion model^ 2 - and the compressibility is 
related to the spinless density of states which has a Id van 
Hove divergence. In contrast the GA yields a compress- 
ibility which tends to zero at the BR point. We remark 
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that the GA compressibility has a jump discontinuity for 
n = 1 and U > U c . In fact, its left and right limits 
are finite, while its value computed in n — 1 is zero. At 
half-filling an antifcrromagnetic (AF) broken symmetry 
TDGA computation instead of the present paramagnetic 
one yields much more accurate results 4& 

In the dilute limit the exact compressibility diverges 
again whereas the TDGA result yields a finite value. This 
disagreement is not surprising since the RPA in general is 
well known to fail at low densities. In this case a particle- 
particle approach, recently implemented on top of the 
GA)£i would be more appropriate. 

Despite the (expected) failure of the paramagnetic 
TDGA at low and half filling at intermediate fillings 
the behavior of the compressibility as a function of U 
is qualitative and to some extent quantitatively repro- 
duced. One should keep in mind again that we are using 
a Fermi liquid approach whereas the real ground state is 
a Luttinger liquid. 



U/U C =0.39GA 
U/U C =0.39QMC 
U/U C =0.78GA 
U/U c =0.78 QMC 




0.6 
0.5 
* 0.4 
0.3 
0.2 



FIG. 3: (Color online) Id-Charge susceptibility (n = 0.6) 
as a function of q/n: comparison between TDGA and QMC 
results.— 

In Fig. [3] we compare the TDGA charge susceptibility 
with QMC results from Hirsch and Scalapino, 19,21 Since 
their data are for n = 0.6 we expect the TDGA to give 
reasonable results. Although our formalism is at T = 
and the QMC study at T ^ 0, the comparison is mean- 
ingful because their results are at T = 0.0690, quite low if 
compared to the electronic energy scales. The QMC sus- 
ceptibilities generally agree with ours within ten-twenty 
percent deviations, and with larger deviations at high q 
for U ~ U c . For large U, in fact, the QMC data exhibit 
the transfer of the peak from q = 2kp to q = 4kp, sig- 
nature of the spin-charge separation of the ld-Luttinger 
liquid, clearly absent in our Id FL. The QMC curves 
present a finite T effect that smoothes the peak. 

The above results indicate that our FL scheme works 
quantitatively rather well away from half-filling (where 
an AF TDGA computation would do a better job) and 
provides reasonable momentum dependencies (but for 
the subtle Luttinger Akp peak shift for large IT). 

In Fig. 21 we show the charge susceptibility Kq as a 
function of U/U c . These results should be taken with a 
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FIG. 4: (Color online) Id-Charge susceptibility as a function 
of U/Uc for n = 0.75. 



pinch of salt due to the explained drawbacks, however 
they illustrate well the general behavior of the TDGA 
that are found in higher dimension where the approach 
performs better^ For small deviations of the density 
from half-filling the compressibility has a minimum close 
to U — U c which is due to the corresponding maximum in 
the (11) element of the interaction kernel (cf. Fig. [I}. For 
large q this minimum becomes too shallow to be clearly 
seen in Fig. [4] At small momenta the charge suscepti- 
bility is close to the compressibility. As the momentum 
approaches q = 2k p — nix the charge susceptibility di- 
verges for small U . However, this divergence is strongly 
suppressed upon increasing U. At small doping n q is fi- 
nite but still shows a shallow minimum close to U = U c . 
This behaviour is due to the proximity of the Mott phase 
which is more clear in higher dimensions^ The essen- 
tial point is that the maximum charge response changes 
from the wave- vector q = 2kp to q = upon increasing 
U. Therefore correlations suppress the nesting induced 
transition to a CDW state in favor of phase separation 
as is discussed in more detailed and higher dimensional 
systems in Ref. [6^. 

It is important to notice that the failures of the TDGA 
found at high filling in our comparison with the exact re- 
sults can be traced back to specific features of the Id 
physics (like, e.g., the spin-charge separation, the equiv- 
alence to spinless fermions at U = oo). Therefore these 
failures should not be attributed to the TDGA per se, 
but to the underlying assumption of a FL ground state. 
This is why our results not only provide qualitative in- 
formations on the ID case, but also shed light on the 
physics of the FL state occurring at higher dimensions. 



B. Holstein Coupling 



The (static) phonon self-energy for the Holstein cou- 
pling S£ o! ee S£ oi (w = 0) = g 2 hol [x q ]u(u = 0) corre- 
sponds to the local charge correlation so that the dis- 
cussion of the previous section directly applies also here. 
Fig. shows Y/l° l for the half-filled system where it is 
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given by 



^1 



2 

9 hoi 



ix° g ] 



ii 



l 



ii 



As is clear from the previous section the primary pur- 
pose of these results is to illustrate the generic behavior 
expected in higher dimensions rather than to comprise 
the physics of ID systems. 

Due to the divergency of [Xq]n &t q = 2kp = ir the 
corresponding singularity in E^°' gets suppressed upon 
increasing U. This suppression persists for all momenta 
which can also be seen from the vertex T q (Fig. [SJd, cf. 
Eq. (|28|1 which quantifies the phonon frequency shift for 
the correlated system as compared to the non-interacting 
case: 



1 



l-4rtx2]n' 



(46) 



In the limit q — > where (— X q )n equals the density 
of states N(Ep) the reduction of T q=0 is therefore de- 
termined by the effective interaction A q which diverges 
for U approaching the Brinkmann-Rice transition (and 
thus T q — > 0). On the other hand, at q = it the lo- 
cal non- interacting charge correlations (Xq=-n-)n display 
a divergence due to nesting and are thus responsible for 
the vanishing of T„^n in this limit. 
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At finite doping (Fig. [6^,) the singularity of T,^° l at 
U = occurs at q — 2kp < it and similar to the half- 
filled system becomes suppressed upon increasing U. As 
discussed in the previous section (cf. Fig. for large 
onsite interaction U the charge susceptibility (and thus 
E^°') acquires a maximum at small momenta so that the 
dominant phonon renormalization is shifted from q = 
2kp to q = 0. On the other hand (Fig. \fijp>) the reduction 
of T q is still most pronounced at the Fermi momenta 
q = 2kp where the bare Lindhard susceptibility (Xq)n 
logarithmically diverges whereas £ g= 2fc F (U > 0) stays 
finite and thus T q= 2k F — 0. Of course this is peculiar to 
the one-dimensional system where one has perfect nesting 
for each carrier density. 
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FIG. 6: (Color online) a) Self-energy E 9 in units of t/ghoi and 
b) ratio between correlated and uncorrelated self-energies T q 
for the Holstein-Hubbard model (particle density n — 0.8) as 
a function of momentum q. 



We now analyze in more detail the mode which be- 
comes soft when the dominating instability shifts from 
q = 2k p to q = for large U and close to half- filling. 
The dressed phonon propagator 



I , 1 , 1 

71/2 t 

q 

FIG. 5: (Color online) a) Self-energy E q in units of t/ghol and 
b) ratio between correlated and uncorrelated self-energies T q 
for the half-filled Holstein-Hubbard model as a function of 
momentum q. 



D < M " TOip (47) 

couples the Holstein phonon Slo with the energy of the 
particle-hole excitations ~ v p q when we use the long 
wavelength limit for the charge susceptibility given in 
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Eq. (O). Then DJu) 



to 



acquires new poles at 



-Q 



(48) 
(49) 



corresponding to a hardening of the Holstein phonon and 
a softening of the effective ('zero sound') particle-hole 
velocity. Thus the instability at q = does not follow 
from a zero in the phonon-type mode but due to the 
fact that the particle-hole excitation acquire a negative 
velocity. Since the poles of the dressed charge propagator 
are identical to those of D q {w) this also corresponds to a 
phase separation instability so that the present approach 
generalizes the analysis of Ref. [l3 for U — > oo Hubbard 
models to finite onsite interactions. 



C. SSH coupling 

For the transitive electron-phonon coupling we have 
seen in Sec. IIIEI that correlations already induce a renor- 
malization of the phonon dispersion Eq. (|39p 



(50) 



due to the elimination of the double occupancy fluctua- 
tions (cf. Eqs. Q33I34P ). Here = 2uj sin(q/2) denotes 
the acoustic branch for U = and the correlation in- 
duced contribution ~ — (Ail q ) 2 is always negative (since 
U q > in Eq. (l3"9")) ). The corresponding softening of f2 9 
is shown in Fig. 
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FIG. 7: (Color online) Acoustic phonon dispersion in the SSH 
model for three different U values at half-filling (g = 0.01, 
u>o = 0.1). Shown is only the influence of the double oc- 
cupancy fluctuations on Q q . The inset displays atomic dis- 
placements with wave- vector q — n/2 (horizontal arrows) and 
the associated modulation of the double occupancy (vertical 
arrows) with the same periodicity. 

The renormalization vanishes for both q — > and 
q — > tt, and is largest for intermediate momenta q ~ tt/2. 



This can be understood from Eq. (|34|) where the first 
term links the displacements to the double occupancy 
fluctuations SD q ~ s'm(q)/U q Q q . Remember that U q is 
the interaction energy of double occupance fluctuations 
(cf. Eq. (|7J) which has a significant momentum depen- 
dence only close to half-filling and large U. Therefore 
the spatial relation between Q q and 6D q is mainly de- 
termined by sin(g) and thus largest at q w tt/2. The 
inset to Fig. [7] depicts the corresponding lattice modula- 
tion (horizontal arrows) which, due to the inreased (de- 
creased) hybridization, favors a modulation of the den- 
sity and double occupancies with the same periodicity 
(vertical arrows). 




FIG. 8: (Color online) Correlation induced correction to the 
phonon dispersion (in units of gt 2 ) from the elimination of 
double occupancy fluctuations. Charge densities are n = 1 
(panel a) and n = 0.8 (panel b) 

The correction (Att q ) 2 to the phonon dispersion in- 
duced by the double occupancy fluctuations is separately 
displayed in Fig. [5) For the half-filled system (top 
panel) the maximum of (Afl q ) 2 shifts to smaller q-values 
upon increasing U due to the more significant momen- 
tum dependence of U q as mentioned above. This is less 
pronounced for the doped system (lower panel) where 
the maximum in the correlation induced correction stays 
close to q — tt/2. Note also that (Afi g ) 2 has a maximum 
as a function of U. This is due to the fact that the SSH 
electron-phonon interaction Eq. §§§ is renormalized by 
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the z-factors which decrease with increasing U so that 
the transitive fluctuations become suppressed. In this 
regard (AQ q ) 2 results from a subtle interplay of kinetic 
and correlation effects. 

We now turn to the influcence of electronic density 
fluctuations on the phonon dispersion which is measured 
in terms of the phonon self-energy S 9 Eq. (|40|) . 
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FIG. 9: (Color online) a) Phonon self-energy E, for the half- 
filled Hubbard-SSH model and different values of U/U c . b) 
Ratio F q = E 9 (t/)/E g ({7 = 0) for the same system. Parame- 
ters: g = 0.01, loo = 0.1. 

Fig. [3^, displays T, q for the half-filled system. In this 
limit the local density fluctuations (originating from the 
Hubbard interaction) are decoupled from the transitive 
ones. Therefore the latter are not screened and the diver- 
gence at q = 7T in case of the SSH coupling is not removed 
upon increasing U in contrast to the Holstcin-Hubbard 
model. As a consequence, the phonon excitations always 
(i.e. for infinitesimally small electron-phonon coupling) 
become instable for q — ir corresponding to the homo- 
geneous dimerized state. Correlations lead to a suppres- 
sion (enhancement) of E 9 for large (small) momenta as 
can be more clearly seen from Fig. ^jp which shows the 
ratio r g = T, q (U) /Y, q (U = 0). In the limit q = ir one 
can show that the ratio is given by T q=7T = Zq, i.e. it is 
completely determined by the hopping renormalization 
factors of the GA. This is consistent with the fact, that 
for the half-filled dimerized system correlations suppress 



the dimerization order paramete r 2 ^ 37 due to the reduc- 
tion of the effective electron-phonon coupling. 

In the limit q — > the self-energy vanishes, however, 
the slope of £ g ^o strongly depends on the correlations 
and leads to the observed increase of r g _>o with increas- 
ing U (Fig. [HJd). The main reason for this enhancement 
comes from the fact that the SSH coupling is a coupling 
to transitive electronic correlations which at half-filling 
and small momenta are decoupled from the local ones. 
In the limit q — ► and half-filling Eq. (|40[) becomes 



^SSH 



s 9 2 l%T(x°) 



22 





\ 2 n U*~ 




H 


U 2 ' 


iq 


^ C _ 


= n 





(51) 
(52) 



and one finds that T, q H is not screened by the strong lo- 
cal charge fluctuations. On the contrary, since (x q )22 ~ 
1/zq it becomes enhanced due to the increase of the 
quasiparticle mass with U. However, for the bare SSH 
coupling \W q \ 2 = q 2 z\ this effect would be overcompen- 
sated resulting in S 9 ~ z 2 . It is due to the TDGA in- 
duced vertex corrections Eq. I]52p that the increase of 
(Xq)22 is even amplified by the concomitant increase of 
\W 2 \ 2 with U. Finally, another (though much weaker) 
factor which leads to the enhancement of S~~q with U 
comes from the dependence of the coupling constant g q 
Eq. (|4ip on the phonon frequencies g q ~ which 
become softened due to the elimination of the double oc- 
cupancy fluctuations (cf. Eq. ([55j)). 

The enhancement of the vertex T q at small momentum 
is a new effect very much in contrast with the result in 
the Holstein case^ 3 - where one always finds T q < 1, i.e. a 
reduction of self-energy corrections with U. 

Fig. HU displays the behavior of S 9 and T q for the 
doped SSH model. Similar to the case of half-filling, E 9 
is reduced upon increasing U for large momenta. How- 
ever, the behavior for small q becomes more subtle as 
can be seen from Fig. [TUb . As a function of U the self- 
energy S g _»o passes through a minimum and for large 
U exceeds again the uncorrelated value (i.e. T q ^o > 1) 
similar to the half-filled case. This behavior results from 
a subtle interplay between local and transitive charge 
fluctuations which are now coupled. For small U the 
screening induced by the local charge fluctuations leads 
to a suppression of (x 9 )22 and also S 9 . Only at larger 
U the vertex corrections for | W 2 \ can overcome this de- 
crease and effectively enhance again the self-energy at 
small momenta. 

The coupling of the local charge density fluctuations 
also contributes to the suppression of the q = 2fcp di- 
vergence in the self-energy. As in case of the Hol- 
stein coupling one therefore finds that T q= 2k F = 
,(U — 0) logarithmically diverges whereas 



since S 9=2fej 
S q= 2fe F (U > 0) stays finite. 

Within the Holstein-Hubbard model we have seen that 
away from half-filling the maximum self-energy shifts 
from the nesting vector q = 2kp to q = when the cor- 
relations become sufficiently strong. As a consequence, 
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FIG. 10: (Color online) a) Phonon self-energy E 9 for the 
doped (n = 0.8) Hubbard-SSH model and different values of 
U/Uc- b) Ratio Y q = E,(f7)/E,(t/ = 0) for the same system. 
Parameters: g — 0.01, u>o — 0.1. 



a large electron-phonon coupling will induce a CDW in- 
stability for small, but a phase separation instability for 
large U , also of course depending on the carrier density. 
Is there a similar scenario for the transitive coupling in 
the SSH model? We determine the instabilities from the 
zero frequency poles of the phonon propagator 



D°(uj = 

DJU = 0) = 9 n/ ; 

which yields the condition 



fin 



(53) 



(54) 



and rt q is the effective phonon dispersion given in Eq. 
(|39|) . The solid line in Fig. QT] marks the lattice in- 
stability at q = 2fci?, i.e. where the system undergoes 
a transition towards a combined CDW and bond-order 
state. This instability is suppressed for large U due to 
the suppression of the 2kF peak in E 9 as shown in Fig. 

nm 

The maximum at n — 0.85 in the instability line is due 
the following. At n — 1 the self-energy can be written 
as S g ({7) = ZqE 9 ([7 = 0) so that upon approaching half- 
filling the self-energy is determined by both the diverging 




FIG. 11: (Color online) Phase diagram for the SSH model 
with parameters g = 0.02, ujq = 0.2. The solid line indi- 
cates the transition to a nesting induced q = 2&f instability 
whereas the dashed-dotted line is the q = instability line. 
The region enclosed by the dashed line corresponds to the 
parameter space where Q. q <Q and the bar at n = 1 indicates 
the localized regime for U > U c . 



E g (/7 = 0) and the Gutzwiller renormalization factor Zq 
which at n = 1 tends to zero for U > U c (in Fig. [TT] the 
solid bar at n — 1 indicates this regime where the charge 
carriers are localized). As a consequence |£ g (£/)| devel- 
ops a maximum as a function of concentration and fixed 
U which is reflected in the maximum of the instability 
line. 

Another instability occurs when the system is stable 
against nesting (i.e. |E g= 2fc F | < ri g =2fe F /2) but the slope 

of | S g ..o I becomes larger than the slope of fi^o- Then 

there exists another solution of the condition Eq. (j5"4")) 
The transition towards this instability occurs at q = 
when both slopes become equal. The corresponding line 
is shown in Fig. [TTj) by the dashed dotted curve. Simi- 
lar to the Holstein- Hubbard model we thus find a q = 
instability for large U which here is confined to a region 
close to half- filling. However, in contrast to the Holstein 
model, where the phase separation is due to an instabil- 
ity of the 'zero sound' particle-hole collective mode v p q 
caused by the coupling to the optical phonon we have in 
the SSH model a coupling between two acoustic modes, 
i.e. Q q v p hq and v p q. For the situation we have ana- 
lyzed in Fig. [Tl]we find always v p h < v p so that the mode 
which becomes unstable has dominantly phonon charac- 
ter. However, since v p in the Mott regime is renormalized 
to very small values one could also imagine a situation 
where v p h > V p in analogy to the Holstein case. 

Finally, we have seen that the TDGA applied to the 
SSH coupling yields an effective phonon dispersion Eq. 

which yields correlation induced softening through 
the elimination of double occupancy fluctuations (cf. Fig. 
|5J|. For large coupling g these frequencies can become 
negative even without the consideration of density fluctu- 
ations. The corresponding regime in Fig. [TT]is enclosed 
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by the dashed line. We find (at least for the present 
model) that this area is always in a parameter regime 
which corresponds to the nesting induced instability and 
therefore never gives rise to a 'real' instability. 



IV. CONCLUSIONS 

We have investigated the renormalization of phonon 
frequencies within the Hubbard-Holstein and Hubbard- 
SSH models based on the TDGA approach. Our consid- 
erations of the Holstein coupling for one-dimensional cor- 
related systems supplements our investigations in higher 
dimensions— and serves as a reference for our computa- 
tions of the SSH coupling. In the latter case we have 
found that correlations influence on the phonon modes 
£l q via two mechanisms. First, the coupling to double 
occupancy fluctuations leads to a softening which has a 
maximum around q = 7r/2, depending on U and dop- 
ing. The second mechanism is the standard screening 
from density fluctuations. However, in this regard our 
TDGA approach goes beyond the standard RPA since it 
incorporates the interaction between phonons and both, 
transitive and induced local density fluctuations. This 
leads to an interesting dependence of the self-energy T, q 
on the local repulsion U since it becomes suppressed for 
large but enhanced for small momenta. 

We have found that also the transitive coupling of the 
Hubbard-SSH model gives rise to an interesting phase 
diagram where correlations can suppress the q = 2fcp 
nesting instability but at the same time are responsible 
for the occurence of a q — instability in the vicinity of 
half-filling. Our calculations therefore indicate that the 
phase separation instability, previously only evidenced 
for Holstein- type couplings, seems to be a generic prop- 
erty of strongly correlated electrons coupled to phonons. 
This is interesting in the context of complex oxides which 
often show nanoscale phase separation. We have re- 
stricted to a short range only model, however the short 
range phase when supplemented with the long range 

Coulomb interaction is well known to lead to mesoscopic 
inhomogeneities^^SiSSieSiZOiZ 1 . 

How do our results apply for higher-dimensional sys- 
tems, especially with regard to the anomalous soft- 
ening of bond-stretching modes in perovskite materi- 
als?— iiL£iL2j2i Consider e.g. the half-breathing mode in 
cuprates which involves the movement of two planar oxy- 
gen ions towards the central Cu ion. The induced change 
of the ionic potential on Cu leads to a Holstein-type cou- 
pling whereas the associated modulation of the Cu-0 
hopping integral gives rise to a coupling of the SSH type. 
Concerning the latter interaction it is interesting that 
the double occupancy induced renormalization in the 2- 
D three-band model would lead to a maximum frequency 
shift at the zone boundary (in contrast to that at q = ir/2 
in the ID SHH model). This kind of interaction there- 
fore induces a downwards dispersion of the half-breathing 
mode which in the lowest approximation just vibrates at 



constant frequency. Obviously, in order to account for the 
doping dependence of the softening one has additionally 
to consider the effect from the density fluctuations enter- 
ing the phonon self-energy. In this regard it would be in- 
teresting to investigate wether our approach can improve 
related Hartree-Fock (HF) calculations within the three- 
band model^ 4 - which give an incorrect (i.e. too small) 
doping dependence of the softening. In fact, since the 
correlation functions in the TDGA incorporate the cor- 
relation induced reduction of the kinetic energy its de- 
pendence on the charge carrier concentration is expected 
to be much more pronounced than in the HF approach. It 
should be noted that calculations of the density response 
for the t J-model also indicate a strong renormalization 
of bond-stretching phonons 7 ! with a larger anomaly oc- 
curing for half-breathing as compared to full-breathing 
modes.— 

Our theory can be easily extended towards ground 
states which break translational symmetry. In this re- 
gard it would be interesting to evaluate the phonon renor- 
malization from striped ground states since there is ex- 
perimental evidence 5 -^ that these textures contribute to 
the anomalous phonon softening at intermediate q val- 
ues in high-temperature superconductors. Since codoped 
LSCO compounds, where static stripe order is unam- 
bigously established, show a rather strong renormal- 
ization it has been argued^ 5 - that the corresponding 
phonon dispersion exhibits a Kohn-type anomaly orig- 
inating from the q = 2kp = tt/2 scattering along the 
half-filled stripes. Since the GA (in contrast to HF) leads 
to half-filled stripes as stable mean-field solutions of the 
Hubbard model*! our approach allows for a test of this 
scenario from a realistic model. 

An interesting issue concerns the question wether the 
correlation induced enhancement of the phonon self- 
energy for small wave-vectors in case of the SSH cou- 
pling also reflects a corresponding enhancement of the 
electron-phonon vertex and eventually superconducting 
correlations. In the present model superconductivity is 
mediated by the transitive lattice fluctuations which, as 
we have shown, are not screened in the same way as 
the local ones and in the long wavelength limit even can 
be underscreened. This scenario thus appears similar in 
spirit to previous works by Capone and coworker s 74 ! 75 
where the enhancement of Cooper-pairing near the Mott 
transition of multiband Hubbard models was investi- 
gated. The interorbital fluctuations in these models are 
only little affected by the correlations (which however 
lead to a decrease of the quasiparticles bandwidth and 
a concomitant increase of the density of states near the 
Fermi level) so that they can effectively enhance super- 
conductivity. An analogous mechanism may also work in 
the correlated SSH model. However, this issue is much 
more involved since the investigation of pair-pair scat- 
tering requires a GA energy functional which is charge- 
rotationally invariant 7 ^ also for the transitive electron- 
phonon coupling Eq. . As a consequence the coupling 
between pair and lattice fluctuations will in general be 
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different from W q eZ ~ p/l given in Eqs. (|37I38|) and will 
considered elsewhere. 
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Appendix 

In the TDGA expansion Eq[l4]we have introduced the 
following abbreviations for the z-factors and its deriva- 
tives: 



= dz ia _ , 

Zia — ^0? o ■ — Z . 
OPiia 

~fy~ = Z+ -> dD~ =Z ° 

d Ziu ft Zi(j '/ Zi(j >> 

d pla ~ ++ ' dpiiadpii-o ~ + ~' dp 2 u _ a ~ — 

d Zia- ft d zip _ a d zip tt 

~dD? = Z °' dp~dD i = Z+Dl dp~~dDi = ^ 

For the half-filled paramagnetic state we have z = 
and z +D = z_ D . 
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